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1 Abstract 

<n : 

We calculate the fraction of B solar neutrinos that arrive at the Earth as a V2 mass eigenstate 
as a function of the neutrino energy. Weighting this fraction with the 8 B neutrino energy spectrum 

00 . 

and the energy dependence of the cross section for the charged current interaction on deuteron with 

a threshold on the kinetic energy of the recoil electrons of 5.5 MeV, we find that the integrated 

weighted fraction of jVs to be 91±2% at the 95% CL. This energy weighting procedure corresponds 

^ \ to the charged current response of the Sudbury Neutrino Observatory (SNO). We have used SNO's 

i ' 

current best fit values for the solar mass squared difference and the mixing angle, obtained by 

<D : 

■ combining the data from all solar neutrino experiments and the reactor data from KamLAND. 



The uncertainty on the v-i fraction comes primarily from the uncertainty on the solar 5m 2 rather 



than from the uncertainty on the solar mixing angle or the Standard Solar Model. Similar results 
for the Super-Kamiokande experiment are also given. We extend this analysis to three neutrinos 
and discuss how to extract the modulus of the Maki-Nakagawa-Sakata mixing matrix element U e 2 
as well as place a lower bound on the electron number density in the solar 8 B neutrino production 
region. 
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I. INTRODUCTION 

Recently the KamLAND and Sudbury Neutrino Observatory (SNO) 2] experiments 
have given a precise determination of the neutrino solar mass squared difference and mixing 
angle responsible for the solar neutrino deficit first observed in the Davis [3J] experiment 
when compared to the theoretical calculations by Bahcall |4|. Subsequently this deficit has 
been observed by many other experiments j^, q|, while the theoretical calculations of the 
neutrino flux based on the Standard Solar Model (SSM) has been significantly improved 3]. 
When all of these results are combined in a two neutrino fit as reported by SNO Q], the 
allowed values for the solar mass squared difference, StUq, and the mixing angle, 8 & , are 
individually (for 1 degree of freedom) restricted to the following range 1 , 

6m% = 8.0+°;^ x l(T 5 eV 2 , 
sin 2 = 0.310 ±0.026, (1) 

at the 68 % confidence level. Maximal mixing, sin 2 ^© = 0.5, has been ruled out at greater 
than 5 a. The solar neutrino data is consistent with v e —>■ and/or v T conversion. The 
precision on StUq comes primarily from the KamLAND experiment l( whereas the precision 
on sin 2 ^ comes primarily from the SNO experiment Q]. 

The physics responsible for the reduction in the solar 8 B electron neutrino flux is the 
Wolfenstein matter effect _M with the electron neutrinos produced above the Mikeyev- 
Smirnov (MS) resonance [lOj. The combination of these two effects in the large mixing 
angle (LMA) region, given by Eq. (fT|). implies that the 8 B solar neutrinos are produced 
and propagate adiabatically to the solar surface, and hence to the earth, as almost a pure 
v 2 mass eigenstate. 2 Since, approximately one third of the v 2 mass eigenstate is z/ e , this 
explains the solar neutrino deficit first reported by Davis. If the 8 B solar neutrinos arriving 
at the Earth were 100% z/ 2 , then the day-time Charged Current (CC) to Neutral Current 
(NC) ratio, CC/NC, measured by SNO would be exactly sin 2 # , the fraction of v e in v 2 in 
the two neutrino analysis. 

Of course, the v 2 mass eigenstate purity of the solar 8 B neutrinos is not 100%, as we will 



1 We use the notation of 8] with the subscript "0" reserved for the two neutrino analysis whereas the 
subscript "12" is reserved for the three neutrino analysis. 

2 Without the matter effect, the fraction of v% 's would be simply sin 2 ^©, i.e. about 31%, and energy 
independent. 
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see later, some fraction arrive as z/^'s and if the electron neutrino has a non-zero component in 
U3 (i.e. non-zero sin 2 # 13 ) then there will be a small fraction arriving as z/ 3 's. For all practical 
solar neutrino experiments, these mass eigenstates can be considered to be incoherent, see 



111 ]. The mass eigenstate purity of the 8 B solar neutrinos is the main subject of this 
paper. In the next section we will summarize the important physics of the MSW-LMA solar 
neutrino solution outlined above and calculate the mass eigenstate purity of 8 B neutrinos as 
a function of the neutrino energy in a two neutrino analysis for both the SNO and Super- 
Kamiokande (SK) experiments. In section 3 we will discuss what happens in a full three 
neutrino analysis. In section 4, as an application of the previous sections, we will discuss 
the possibility of extracting information about the solar interior independently from the 
standard solar model. Finally, in section 5, we present our summary and conclusions. 

II. TWO NEUTRINO ANALYSIS: 
A. 8 B V2 Fraction 

In the two neutrino analysis, let fi(E u ) and / 2 (-£v) be the fraction of 8 B solar neutrinos 
of energy E v which exit the Sun and thus arrive at the Earth's surface as either a v\ or a 
v 2 mass eigenstate, respectively. Following the analytical studies of Ref. these fractions 
are given by 

h{E v ) = (cos 2 ^-P 1 .cos2^) 8B , (2) 
f 2 (E u ) = (sm 2 e% + P x co826%)* B , (3) 

where 6q is the mixing angle defined at the v e production point, P x is the probability of the 
neutrino to jump from one mass eigenstate to the other during the MS-resonance crossing, 
and the sum is constrained to be 1, fi + f% = 1. The average (• • • )sb is over the electron 
density of the 8 B v e production region in the center of the Sun predicted by the Standard 
Solar Model Q]. The mixing angle, 9^, and the mass difference squared, 5m 2 N) at the 
production point are 



sin 2 ^ = - [ 1+ (A - 5 nip cos 29 Q ) _ , 

2[ y/fSmfe cos 26 e - A) 2 + (5m| sin 2# ) ? ' ' 



5m 2 N = J(5m 2 Q cos 26 & - A) 2 + {5m% sin 29 Q ) 2 (5) 
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FIG. 1: The solid and dashed (blue) lines are the 90, 65, 35 and 10% iso-contours of the fraction 
of the solar 8 B neutrinos that are zVs in the drriQ and sin 2 9 & plane. The current best fit value, 
indicated by the open circle with the cross, is close to the 90% contour. The iso-contour for 
an electron neutrino survival probability, P ee , equal to 35% is the dot-dashed (red) "triangle" 
formed by the 65% V2 purity contour for small sin 2 6q and a vertical line in the pure V2 region at 
sin 2 9q = 0.35. Except at the top and bottom right hand corners of this triangle the v% purity is 
either 65% or 100%. 



where 



A = 2V2G F (Y e p/M n )E u = 1.53 x 10" 4 eV 2 f YeP _f" \ 

\kg.cm MeV / 



■ . (6) 

,kg.cm- 3 MeV J W 

is the matter potential, E v is the neutrino energy, Gf is the Fermi constant, Y e is the electron 
fraction (the number of electron per nucleon), M n is the nucleon mass and p is the matter 
density. The combination Y e p/M n is just the number density of electrons. 
Fig. shows, for a wide range of SttIq and sin 2 9@, the iso-contours of 



h = (f2(E v )) E , 



(7) 



where (•••)£ is the average over the B neutrino energy spectrum [14( convoluted with 



□ 



the energy dependence of the CC interaction u e + d^p + p + e cross section LL5II at 
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SNO with the threshold on the recoil electron's kinetic energy of 5.5 MeV. Here we use 
sin 2 #0 as the metric for the mixing angle as it is the fraction of z/ e 's in the vacuum v 2 mass 
eigenstate. In this work, we mainly focus on SNO rather than SK since the former is the 
unique solar neutrino experiment which can measure the total active 8 B neutrino flux as well 
as 8 B electron neutrino flux, independently from the SSM prediction and other experiments. 
However, we give a brief discussion on SK later in this section. 

'. n the LMA region the propagation of the neutrino inside the Sun is highly adiabatic 
16^ . i.e. P x ~ 0, therefore, 
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h{E v ) = l-h{E v ) = (sinXK. (8) 

Due to the fact that 8 B neutrinos are produced in a region where the density is significantly 
higher (about a factor of four) than that of the MS-resonance value, the average (f2(E v ))E 
is close to 90% for the current solar best fit values of the mixing parameters from the recent 
KamLAND plus SNO analysis |2j. Since sin 2 6q — > 1 when A/8m 2 & — > oo (see Eq. (0}), we 
can see that at the high energy end of the 8 B neutrinos (sin 2 9q)»b must be close to 1. 

We can check our result using the analysis of SNO with a simple back of the envelope 
calculation. In terms of the fraction of v% and v 2 the day-time CC/NC of SNO, which is 
equal to the day-time average v e survival probability, (P ee ), is given by 

CC 
NC 

where fi and f 2 are understood to be the V\ and v 2 fractions, respectively, averaged over 
the 8 B neutrino energy weighted with the CC cross section, as mentioned before. Using the 
central values reported by SNO 3 , 

CC 



= (P ee ) =f lC os 2 e G + f 2 sm 2 e Q , (9) 

day 



= 0.347 ±0.038, (10) 

day 



NC 

which was obtained from Table XXVI of Ref. |2j, and the current best fit value of the mixing 
angle, we find f 2 = (1 — fi) ~ 90%, as expected. Due to the correlations in the uncertainties 



3 For the sake of simplicity and transparentness of the discussion, we have avoided the Earth matter effect 
which causes the so called regeneration of v e during night, by simply restricting our analysis to the day 
time neutrino flux throughout this paper. We note that due to the large error, the observed night- 
day asymmetry at SNO is consistent with any value from -8 to 5% |2J whereas the expected night-day 
asymmetry, 2(N-D)/(N+D), is about 2.2-3.5% for the current allowed solar mixing parameters |21|. Thus 
the difference between the day and the day plus night average CC/NC is less than 2% and much smaller 
than SNO's 10% measurement uncertainty on CC/NC. 
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FIG. 2: SNO's Day-time CC/NC ratio in the Suiq versus sin 2 9q plane. At small values of 5m|, 
the Day-time CC/NC ratio equals sin 2 6q. The current allowed region at 68 and 95% CL from the 
combined fit of KamLAND and solar neutrino data 2] are also shown by the shaded areas with 
the best fit indicated by the star. 

between the CC/NC ratio and sin 2 9 & we are unable to estimate the uncertainty on f 2 here. 
Note, that if the fraction of v 2 were 100%, then ^ = sin 2 9 Q . 
Alternatively, we can rewrite Eq. as 4 



Thus how much CC/NC differs from sin 2 9 & is determined by how much f 2 differs from 100%, 
i.e. the size of f\. In Fig. [21 we have plotted the contours of the day-time CC/NC ratio in 
the sin 2 9 & versus SitIq plane for the LMA region. Clearly, at smaller values of StUq the day 
time CC/NC tracks sin 2 9 & whereas at larger values an appreciable difference appears. This 
difference is caused by a decrease (increase) in the fraction that is v 2 {v\) as 5m 2 & gets larger. 
Hence if we know the i>\ or v 2 fraction we can easily calculate sin 2 9 Q from Eq. (fTTj) using a 



4 The relationship between day-time and 6*© (= arcsiny — /i)/(l — 2/i)J or tan 2 6> 

(= — /i)/(l — fi — §§)) is not as transparent as sin 2 8 Q . 




(11) 
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measured value of the day-time CC/NC ratio. 

A similar analysis can also be performed using the event rate of the elastic scattering 
(ES) at SK and/or at SNO. In fact, ES is related to the V\ and v 2 fractions, as follows, 

H = A(cos 2 6 Q + r sin 2 6 Q ) + / 2 (sin 2 6 & + r cos 2 6 & ) (12) 

where r = (a v a, T e) / {cv e e} ~ 0.155 is the ratio of the ES cross sections for u^ T and u e jl^ . 
averaged over the observed neutrino spectrum. Note that we are normalizing the ES event 
rate to that of SNO NC such that Eq. (|T2*j) is valid independent of the SSM prediction of 
the 8 B neutrino flux. 

In general, in the presence of neutrino flavor transitions, the fraction of V\ and v 2 are not 
the same for ES and CC because the energy dependence of the cross sections are different. 
However, in Ref. it was suggested that if we set analysis threshold energies for SK and 
SNO appropriately as Tsno = 0.995 Tsk — 1-71 (MeV), where Tsno and Tsk are the kinetic 
energy threshold of the resulting electron, the energy response of these detectors become 



practically identical 



In]. 



Thus, using such a set of thresholds, even if there is a spectral 
distortion in the recoil electron energy spectrum, to a good approximation, SK/SNO ES and 
SNO CC are related as follows, 

ES CC / CC\ 

NC = NC + T-NCJ' (13) 

and all the results we obtained for SNO in this paper are equally valid for ES at SK and/or 
at SNO provided the energy thresholds are set appropriately 5 . 

In Fig. [Hl^a) we show the v 2 fraction, f 2 (E u ), versus E v . The rapid decrease in the v 2 
fraction below E v ~ 8 MeV is responsible for the expected spectral distortion at energies 
near threshold in both SNO (see Fig. 36 of Ref. |2j) and SK (see Fig. 51 of the last Ref. 
in P]). For a neutrino energy near 10 MeV, the SNO sweet spot, the 90% CL variation in 
5m 2 & changes f 2 (E u ) more than the 90% CL variation in sin 2 # Q . Whereas in Fig. E£b) we 
give the fraction of u 2 s above a given energy both unweighted and weighted by the energy 
dependence of the CC interaction and ES cross sections. Note, that above a neutrino energy 
of 7.5 MeV there is little difference between the weighted and unweighted integrated v 2 
fraction. Furthermore, in Fig.|3fc), we show the fraction of v 2 § above a given kinetic energy 



5 In fact this suggest an alternative to looking for a spectral distortion to test MSW, compare ES to (1-r) 
CC + rNC for a variety of kinetic energy thresholds. 
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FIG. 3: (a) The fraction of V2, f2(E u ), as a function of the neutrino energy. The solid (black) curve 



is obtained using the central values for 5m,Q = 8.0 x 10 5 eV 2 and sin 2 6q = 0.31 whereas the blue 



dashed (red dotted) lines are the 90% CL range varying StUq (sin 2 6q) but holding sin 2 9 & (5m@) 
fixed at the central value, Eq. (fT|l. (b) The integrated fraction of 8 B neutrinos which are iVs above 
an energy, E u , dashed (red) curve. Whereas, the solid black and blue curves are weighted by the 
energy dependence of the charge current (CC) cross section |l5j ] and the elastic scattering (ES) 
cross section jl8|. respectively, (c) The integrated fraction of 8 B neutrinos as a function of the 
threshold kinetic energy of the recoil electrons for CC (SNO) and ES (SK or SNO) reactions. 



for the recoil electron for both CC (SNO) and ES (SK or SNO) reactions. We observe that 
for the same threshold, f'2 for ES is always smaller than that for CC. This is expected since 
unweighted f2 is a increasing function of E v and CC cross section increase more rapidly with 
energy than that of ES cross section. Hereafter, unless otherwise stated, we focus on the 
SNO CC reaction, as the results for ES reaction are qualitatively similar and the thresholds 
can be adjusted to give identical results for all practical purposes. 

In Fig. E] we give the breakdown into z/ x and z/ 2 for the raw 8 B spectrum as well as the 
spectrum weighted by the energy dependence of the CC interaction using a threshold of 5.5 
MeV for the kinetic energy of the recoil electrons. Here we have used the current best fit 
values for 5m 2 Q and sin 2 9q . 

How does the fraction of z/ 2 vary if we allow 5m 2 Q and sin 2 6 & to deviate from their best 
fit values? In Fig. a) we show the contours of the fraction of v-i in the Sm 2 ^ versus sin 2 9 & 
plane where we have weighted the spectrum by the energy dependence of the CC interaction 
cross-section, and we have used a threshold on the kinetic energy of the recoil electrons of 5.5 
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FIG. 4: The normalized 8 B energy spectrum broken into the v\ and components. The left hand 
curves (black and white) are unweighted whereas the right hand curves (blue and red) are weighted 
by the energy dependence of the CC cross section Q| with a threshold of 5.5 MeV for the recoil 
electron's kinetic energy. 

MeV. This energy dependence mimics the energy dependence of the SNO detector. Because 
of the strong correlation between sin 2 Q and the day-time CC/NC ratio we also give the 
contours of the fraction of v 2 in the 5m 2 & versus day-time CC/NC plane in Fig. Efb). Thus 
the 8 B energy weighted average fraction of zVs observed by SNO is 



/ 2 = 91 ± 2% at the 95% CL. 



(14) 



This is the two neutrino answer to the question posed in the title of this paper. We note, 
however, that as we showed in Fig. |Hfc) the value of f 2 is a function of the threshold energy 
and also depends on the experiment. We estimate that for SK with the current 4.5 MeV 
threshold for the kinetic energy of the recoil electrons, that 



f 2 = 88 ± 2% at the 95% CL. 



(15) 



The uncertainty is dominated by the uncertainty in 5m? & /A. However, the uncertainty on 
SttIq is approximately 5% from the KamLAND data whereas the uncertainty on the matter 
potential, A, in the region of 8 B production of the Standard Solar Model is 1-2%, see 
Hence, the uncertainty on 5m 2 & dominates. 
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B: v 2 fraction (%) 
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FIG. 5: (a) The v 2 fraction (%) in the Sthq versus sin 2 9 & plane. As in Fig. [21 the current allowed 
region is also shown, (b) The v 2 fraction (%) in the SuIq versus the Day-time CC/NC ratio of 
SNO plane. We have excluded a region in the top left hand corner of this plot which corresponds 
to sin 2 6> < 0.1. The current allowed range is indicated by the cross. 
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For the current allowed values for 5m 2 & and sin 2 ^, the ratio 



SttIq sin 28 G 



'0j 



3 
4' 



(16) 



A( 8 B) -8m% cos 2^ 

where A( 8 B) is obtained using a typical number density of electrons at 8 B neutrino produc- 
tion (Y e p pa 90 g.cm~ 3 ) and the typical energy of the observed 8 B neutrinos (~ 10 MeV). 

For the best fit central values of 8m 2 Q and sin 2 # , given by Eq. Q , let us define an effective 
matter potential for the 8 B neutrinos, A 8 ^*, such that the left hand side of Eq.(jlJ) equals 
our best fit value for the fraction that is v 2 . Thus, 

A*£ f = 5m 2 & sin 26 Q cot 26> + 



2/ 2 -l 



2v7 2 (i-/ 2 ) 



(17) 



1.36 x 10" 4 eV 2 
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for f 2 = 0.910. This A e j f corresponds to a Y e pE v = 0.892 kg cm MeV, the effective 
mixing ang le, e%\ eff = 73° and the effective 5m 2 N \ eff = 13.6 x 10~ 5 eV 2 . 

We can then use this ^4*// ^° perform a Taylor series expansion about the best fit point 
as follows 

A = <- 2 o.B^ + |U+ (a with {s |- us, 

10 125 4 {A e ff — (wig cos 20 Q ) 

This simple expression reproduces the values of to high precision throughout the 95% 
allowed region of the KamLAND and the solar neutrino experiments given in Fig. EJa). In 
this sense our ^4*// * s the effective matter potential for the 8 B neutrinos. An expansion in 
5m G /A around its typical value of 0.6 could also be used but the coefficients are ever more 
complex trigonometric functions of 9®, whereas with our £ expansion the coefficients are 
small rational numbers. 



B. 'Be and pp neutrinos 

For 7 Be and pp neutrinos the fractions of v\ and v<i are much closer to the vacuum values 
of cos 2 #0 and sin 2 & respectively, as they are produced well below (more than a factor 
of two) the MS-resonance in the Sun, and an expansion in A/8m & is the natural one. In 
the third Ref. in Q], the electron neutrino survival probability was obtained by a similar 
expansion around the average of the matter potential. Using this expansion, we find that 

/2=1 _ /l=sin * c=sinX + i sill2Me (_^) +0 (_y 2 (19) 

with A 7 e f f = 1.1 x 10" 5 eV 2 and Af f} = 0.31 x 10" 5 eV 2 , (20) 

where the averaged value of the energy (weighted by the cross section) as well as the electron 
densities used are, respectively, (E u ) pp = 0.33 MeV and (Y e p) pp = 62 g/cm 3 for pp, and 
(E u )r Be = 0.86 MeV and (V e p)7 Bc = 81 g/cm 3 for 7 Be. Thus / 2 ( 7 Be) = 37 ± 4(7)% and 
/2(pp) = 33 ± 4(7)% at 68 (95) % CL where the uncertainty here is dominated by our 
knowledge of sin 2 6 & . 



C. Two Neutrino Summary 

In Fig. El we give the neutrino mass spectrum, the value of fraction of zVs (sin 2 ^) and 
the fractional flux as function of the electron number density times neutrino energy, Y e pE u , 
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Y e pE v (kg cni Me V) 

FIG. 6: The Mass spectrum (top panel), the fraction of iVs produced, sin 2 OS , (middle panel) and 
the fractional flux (bottom panel) versus the product of the electron fraction, Y e , the matter density, 
p, and the neutrino energy, E u , for the best fit values SvOq = 8.0 x 10~ 5 eV 2 and sin 2 9 Q = 0.310. 
The vertical dashed lines give the value of Y e pE v which reproduces the average v<i fractions, 91, 37 
and 33% for 8 B, 7 Be and pp respectively. This value of Y e pE u = 0.89 kg cm -3 MeV, for the 8 B 
neutrinos, gives a production mixing angle equal to 73° and a production $m 2 N = 14 x 10~ 5 eV 2 . 
Y e pE v = 1 kg cm -3 MeV corresponds, in terms of the matter potential, to 15.3 x 10~ 5 eV 2 , see 
Eq. ©. 



which is proportional to the matter potential, for the 8 B, Be and pp neutrinos using the 
best fit values of dm^ and sin 2 9 Q in Eq. The 8 B energy spectrum has been weighted 
by the energy dependence of the CC interaction of SNO with a 5.5 MeV threshold on the 
kinetic energy of the recoil electrons whereas the pp energy spectrum has been weighted 
by the energy dependence of the charged current interaction on Gallium with a 0.24 MeV 
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threshold. The vertical dashed lines gives the value of Y e pE v which reproduces the average 
v 2 fraction using the simple expression in Eq. (j3J) and are useful for the approximations given 
in Eqs. (jTHj) and lfI3). 

The energy weighted z/ 2 fractions for 8 B, 7 Be and pp neutrinos using a two neutrino 
analysis, at the 95% CL, are 

/ 2 ( 8 B) = 91 ± 2%, (21) 
/ 2 ( 7 Be) = 37 ±4%, (22) 
/ 2 (pp) = 33 ±4%, (23) 

where the uncertainties for 7 Be and pp are dominated by the uncertainty on sin 2 & whereas 
for 8 B the uncertainty is dominated by the uncertainty on 5m 2 & . The v\ fractions, fi, are 



simply 1 - /; 



2- 



III. THREE NEUTRINO ANALYSIS 

For the three neutrino analysis we first must discuss the size of the component of u 3 
which is z/ e , i.e. the size of sin 2 #13. This mixing angle determines the size of the effects on v e 
associated with the atmospheric mass squared difference. The best constraint on #13 comes 
from the CHOOZ reactor experiment [2Cj which gives a limit on sin 2 #13, as 

< sin 2 9 13 < 0.04, (24) 

at the 90 % CL for Sm^ = 2.5 x 10~ 3 eV 2 . This constraint depends on the precise value of 
Sm^i with a stronger (weaker) constraint at higher (lower) allowed values of Sm^. 

So far the inclusion of genuine three flavor effects has not been important because these 
effects are controlled by the two small parameters 

— |i « 0.03 and/or sin 2 13 < 0.04. (25) 
dm 32 

However as the accuracy of the neutrino data 



account genuine three flavor effects. See 
solar neutrinos. 



21, 



imp roves it will become inevitable to take into 



22 



for recent studies on the impact of #13 on 



Suppose that Double CHOOZ [23j, T2K [2J or NOi/A |25j or some other experiment 
measures a non-zero value for sin 2 #13. What effect does this have on the previous analysis? 
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How does this change our knowledge of the solar parameters and the relationship between 
solar mixing angle and the fraction of v<p. 

Our knowledge of the solar 5m 2 comes primarily from the KamLAND experiment where 
the effects of the atmospheric Sm 2 are averaged over many oscillations, thus to high accuracy 

fi m 21 = ^ m 0; (26) 

i.e. the solar Sm 2 remains unaffected. Remember, we are using the notation 5m 2 i and 
sin 2 9i2 for the three neutrino analysis to distinguish it from 5m 2 & and sin 2 used in the 
two neutrino analysis. 



A. 8 B 3 Neutrino Analysis 

For the mixing angle sin 2 9 12 the situations is more complicated in the three neutrino 
analysis. The 8 B electron neutrino survival probability measured by SNO's day-time CC/NC 
ratio can be written as 

CC 

— = T\ cos 2 6>i3 cos 2 0i2 +T2 cos 2 0i3 sin 2 9 12 + F 3 sin 2 9 13 , (27) 

where T\, T2 and T 3 are the fraction of V\, v 2 and u 3 respectively, satisfying JF X +JF 2 +JF 3 = 1. 
The U3 fraction is given by 

T z = (l ± j^^j sin 2 0i3 « sin 2 9 13 , (28) 

where +(-) sign refers to the normal, 8m\ x > (inverted, Sm 31 < 0) mass hierarchy. The 
small correction factor , 2 \ , ~ 10% comes from matter effects associated with atmospheric 

l" m 31 1 

Sm 2 in the center of the Sun. We will ignore this correction since it is small and currently 
the sign is unknown. Hence, JF X + JF 2 = 1 — JF 3 = cos 2 0i3. 

With this approximation the v\ and z/ 2 fractions can be written as 

;Fi = cos 2 0i3 (cos 2 0^)s B and T 2 = cos 2 9 13 (sin 2 0^) 8B . (29) 

where the average (• • • )sb is over the solar production region and the energy of the observed 
neutrinos, sin 2 9^ 2 is given by Eq. (J3J) with the replacements sin 2 — > sin 2 i2 and A 

A COS 2 013 

In going from the two neutrino analysis to the three neutrino analysis the quantity that 
must remain unchanged is the value of the electron neutrino survival probability, i.e. the 
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CC/NC ratio. This implies that we must adjust the value of sin 2 9 12 and hence the fractions 
of V\ and v 2 so that the CC/NC ratio remains constant. We have performed this procedure 
numerically and report the result as a Taylor series expansion in the fraction of zVs about 
sin 2 9i3 = 0. If we write 

^(sin 2 0i 3 ) = ^i(O) + a sin 2 9 13 + O (sin 4 13 ), (30) 



then .Fi(O) = /1, and a = , 

a sin 



13 



In FigEfa) we have plotted the contours of a = ds ^\ 
Near the best values this total derivative is close to zero, i.e. 



(31) 

in the 5m 2 & versus sin 2 9 & plane. 



d sin 2 



13 



= 0.00^ 2 (32) 

sin 2 0i3=O 



at the 68% CL. As sin 2 0i3 grows above zero, the size of T\ is influenced by a number of 
effects; the first is the factor of cos 2 i3 in Eq. (|2*U|) which reduces the second is the 
matter potential A which is reduced to A cos 2 0i 3 raising the fraction T\ and third is the 
value of sin 2 0i2 which changes to hold the CC/NC ratio fixed. By coincidence the sum 
of these effects approximately cancel at the current best fit values and the fraction of v\ 
remains approximately unchanged as sin 2 0i 3 gets larger. This implies that the fraction of 
Vi is reduced by ~ sin 2 i3 since the sum of T\ + T% is simply cos 2 9is, thus 

F x » f x = 0.09 =F 0.02, (33) 
F 2 = f 2 - sin 2 13 ss 0.91 ± 0.02 - sin 2 13 , (34) 
^3 = sin 2 13 . (35) 

Remember fa and Ti are the fractions of the i-th mass eigenstate in the two and three 
neutrino analysis, respectively. The uncertainty comes primarily from the uncertainty in 
SuIq measured by KamLAND. 

As a use of these fractions one can for example evaluate the MNS matrix element, \U e 2\ 2 = 
cos 2 0i3 sin 2 0i2, by rewriting Eq. ()27|1 as 

\tt |2 2 Q -2a (i§ _ cos2 ^13-^1) / Qft x 

CLo = cos 0i3 sin 0i2 = , - ^r^, 36 

1 ei| u (cos 2 13 - 2T X ) ' v ; 

where terms of O(sin 4 13 ) have been dropped. Performing a Taylor series expansion about 
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o^dF/dsiire^ l (%) 



MIU e2 | 2 /dsin 2 e 13 | (%) 





sin 2 



sin 2 9 r 



FIG. 7: Iso-contours of the derivatives of T\ (a) and |?7 e2 | 2 (b) with respect sin 2 #13 evaluated at 
sin 2 #13 = in the <5mg versus sin 2 6 & plane. The contours are labeled as in per cent. The 68 and 
95 % CL allowed regions are also indicated. 



sin 2 6*13 = 0, we find 



\U, 



c2\ 



sin 2 + 9 13 + (sin* 9 



13; 



with (3 



d\U, 



c2 



d sin 2 #13 



(ft - a) + (l + 2o;)sm 2 







(l-2/i) 



For the current allowed region of the solar parameters, this implies that 



|[/ e2 | 2 « sin 2 ^ + (0.53+H4)sin 2 #i 3 



(37) 
(38) 

(39) 



at the 68% CL, i.e. the three neutrino |£4 2 | 2 i s approximately equal to the sin 2 #q B using a two 
neutrino analysis of only the 8 B electron neutrino survival probability using the KamLAND's 
5rriQ constraint plus 53% of \U e 3\ 2 determined, say, by a CHOOZ-like reactor experiment, 
see Fig. E^b). 

If a similar analysis is performed for the three neutrino sine squared solar mixing angle 
sin 2 #12, the total derivative with respect to sin 2 613 is simply {(3 + sin 2 9 & ). For tan 2 612 the 
total derivative is (/3 + sin 2 9 Q ) / cos 4 9 Q . Alternatively we can turn this discussion inside out 
and write the 8 B effective two component sin 2 Q in terms of three component quantities as 

sin 2 #q B = sin 2 12 - (J3 + sin 2 12 ) sin 2 6 13 . (40) 
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For KamLAND, the equivalent relationship is 

^ C » = sin ^ 12 -(|^) sin ^ 13 . (4!) 

For the current best fit values (/3 + sin 2 6 l i2) ~ 0.90 is close to sin 2 2#i2 /2 cos 2#i2 ~ 1.1, i.e. 
in a two component analysis the difference between the solar 8 B and KamLAND sin 2 # 's is 
approximately 0.2 sin 2 #13. 



B. 7 Be and pp 3 Neutrino Analysis 

Performing a similar 3 neutrino analysis for the pp (or 7 Be) neutrinos we find that the 
fraction of neutrino mass eigenstates is 

T x « cos 2 6 & - I sin 2 26 Q (-^ + sin 2 9 l3 =fi + 0.82 sin 2 6 13 , (42) 

L \^07TIq J COS 2uq 

T 2 « sin 2 9 Q + I sin 2 26 Q ( -^-] - sin 2 13 = f 2 - 1.8 sin 2 13 , (43) 



2 \ 5m| / cos 26 & 

^3 ~ sin 2 fl 13 , (44) 

where the sin 2 9 Q here is determined from the pp (or 7 Be) neutrinos. Terms of order 
O [A/SttIqY , 0(sin 4 13 ) and O (sin 2 6i 3 A/ Sm^j have been dropped here. The two neu- 
trino fractions fx and f'2 are given in Eq. ()19|) . 

Again we can use these fractions to determine the \U e 2\ 2 element of the MNS matrix 

\U e2 \ 2 = sm 2 6 Q - j sin 2 # 13 « sin 2 # - 1.8 sin 2 fl 13 . (45) 



\ cos 2^0 / 

Comparing this equation with Eq. (J3~9"j) appears to be in contradiction but this is not so 
since if sin 2 #13 7^ then the two component analysis of the 8 B and pp (or 7 Be) neutrinos 
will lead to different values of sin 2 # , in fact 

sin 2 6™ - sin 2 fl B « 2.3 sin 2 6 13 . (46) 



This difference has been extensively exploited in Ref. 22j to determine sin 2 # 13 using only 
solar neutrino experiments. Their sin 2 (9 12 versus sin 2 # 13 figures, e.g. Fig. 6, demonstrates 
this point in a clear and useful fashion. Also, the numerical values of our derivatives of \U e 2\ 2 
are consistent with the inverse of the slopes of their Fig. 6. 

Eqs. ()39|) and (|4"5Jl also imply that the uncertainty in the determination of \U e 2\ 2 from 
the current unknown value of sin 2 #i 3 is smaller for the analysis of 8 B neutrinos than pp or 
7 Be neutrinos. Of course the current uncertainty on the two neutrino sin 2 6 & dominates. 
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IV. PROBING THE SOLAR INTERIOR BY 8 B NEUTRINOS 



In this section, as an application of our analysis, we will invert the discussions found in 
Ref. where the validity of the MSW physics has been tested assuming the standard solar 
model (SSM) prediction of the electron number density as well as 8 B neutrino production 
region. Here, we will discuss what can be said about these quantities, assuming the validity of 
the MSW effect in the LMA region. While there is no strong reason to doubt the correctness 
of the SSM, which is in good agreement also with the helioseismological data it is 
nevertheless interesting if we can test it independently. 

Since the propagation of 8 B neutrinos, in the Sun, is highly adiabatic in the LMA region, 
the fraction of z/ 2 , and consequently, the SNO CC/NC ratio is determined only by the 
effective value of the matter potential, defined in Section 11(A). This implies that if 

we can measure sin 2 9 & using an experiment independent of the 8 B solar neutrinos, then 
from the measured value of SNO's CC/NC ratio we can determine the value of ^4 e //- Note, 
that we can not extract information on the electron number density distribution or the 8 B 
neutrino production distribution, separately, but only on ^4g// which is a single characteristic 
of the convolution of these two distributions. 

For the two flavor neutrino analysis, if we rewrite the definition of the effective matter 
potential given by Eq. ffTTj) using the relationship between f 2 and SNO's CC/NC ratio, 
Eq.flHJ), we obtain 

1 - 2-F c 



Ajjf = 5m 2 & sin 26» 



cot 26r, + ~ NC 



2 A /(cos^ -gg)(gg -sin 2 # ) 



(47) 



This expression allows us to obtain a value of ^4*// fr° m (sin 2 9@, 5m 2 & ) measured independent 
of 8 B neutrinos and SNO's 8 B neutrino CC/NC ratio. We can convert this into an effective 
value of the electron number density, Y e p in the solar 8 B production region, as follows 

s B _ M n A e ] f 

Yep Lff = 27^; <^>v (48) 

where {E v )s B = 10.5 MeV is the CC cross section weighted average energy of neutrinos 
observed by SNO. For a given solar model, the value of Y e p f^j- can be calculated for any 
value of sin 2 and 5m 2 & . The SSM prediction is that Y e p = 85 g cm -3 at the current 
best fit point 6 . As a comparison the mean value of Y e p over the 8 B production region is 90 

6 Because of the way we have defined A^^, our Y e p |^ has a weak dependence on sin 2 8 & and Sttlq but 
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FIG. 8: The iso-contours of Y e p in the sin 2 9 & — CC/NCjday plane. The line labeled SSM is 
the Standard Solar Model prediction for Y e p f^j- (~ 85 g/cm 3 ). The range of observed values of 
CC/NC are indicated by the shaded horizontal bands. The KamLAND experiment places a lower 



D. 



bound on sin 0© independent of solar neutrinos at 0.17, see []J]. The vertical band indicate the 
uncertainty which could be expected by future reactor experiments 



g cm -3 . The reason that Y e p f^j- is below the mean value is because values of Y e p below the 
mean pull down the v 2 fraction more than values above the mean raise the v 2 fraction. 

We show in Fig. |HJ the iso-contours of Y e p |^ in the sin 2 9 Q — CC/NC|d ay plane, for the 
current best fitted value of 8m@. The observed range of SNO's CC/NC are shown by the 
horizontal lines 7 . From this plot, we can derive the lower bound on Y e p \ s e ^ which is 40 
g/cm 3 for any value of 9 Q at 95 % CL. Future reactor neutrino oscillation experiments j^l 
can perform a 2-3% measurement of sin 2 9 Q . The 68% range of sin 2 9 & is indicated by vertical 
lines in this figure. However, such precision on sin 2 9 & will not reduce the allowed values for 
Y e p unless the error on the measured value of CC/NC is reduced. 

A three neutrino analysis is needed if U e s ^ and this can be performed using Eq. ()47|) 



this variation is less than 2% over the 95% CL allowed region. 
7 Another horizontal band could be included by combining the Super-Kamiokande Electron Scattering 
measurement with the SNO Neutral Current measurement. However, since the uncertainty on the NC 
measurement dominates this would produce a similar sized band. 
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with the following replacements, 

CC CC 1 

6 Q ^6 12 , 5m 2 Q -f 5m 2 2l / cos 2 6 13 and — -> — cQs4g • (49) 

A weak upper bound could be derived using a precision measurement of the 7 Be and/or pp 
electron neutrino survival probability in a similar fashion. As A e ff gets larger, the fraction of 
v-i gets larger, see Eq.(|19|). and hence the electron neutrino survival probability gets smaller 
for fixed values of the mixing parameters. The upper bound arises when this survival 
probability is below the measured survival probability at some confidence level, assuming 
that the mixing parameters have been determined independent of these solar neutrinos. 



V. SUMMARY AND CONCLUSIONS 

We have performed an extensive analysis of the mass eigenstate fractions of 8 B solar 
neutrinos using only two mass eigenstates (sin 2 8 13 = 0) and with three mass eigenstates 
(sin 2 #i3 7^ 0). In the two neutrino analysis the z/2-fraction is 91 ± 2%. The remaining 9 =F 
2% is, of course, in the v\ mass eigenstate. With these fractions in hand, which are primarily 
determined by the solar 5m 2 measured by the KamLAND experiment, the sine squared of 
the solar mixing angle is simply related to CC/NC ratio measured by the SNO experiment. 
For completeness the mass eigenstate fractions for 7 Be and pp are also given. 

Allowing for small but non-zero sin 2 #13, in a full three neutrino analysis, we found very 
little change in the fraction of Ui's. This implies, since the v 3 fraction is sin 2 13 , that the v% 
fraction is reduced by sin 2 #13. That is, the ^-fraction is 

91 ± 2 - 100 sin 2 9 13 % at the 95% CL. (50) 

Since the CHOOZ experiment constrains the value of sin 2 #13 < 0.04 at the 90% CL this 
places a lower bound on the v 2 fraction of 8 B solar neutrinos in the mid-eighty percent range 
making the 8 B solar neutrinos the purest mass eigenstate neutrino beam known so far, and 
it is a v-i beam! 

As an example of the use of these mass eigenstate fractions, we have shown that for the 
8 B neutrinos observed by the SNO experiment, the [/^-element of the MNS matrix is given 
by 

\U e2 \ 2 « sin 2 6% B + (0.53i°;°4) sin 2 6 l3 . (51) 
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Where sin 9* is the sine squared of the solar mixing angle determined by using a two 
neutrino analysis of the 8 B neutrinos plus KamLAND. An analysis for this sin 2 9^ obtained 
from the SK, SNO and KamLAND data Q gives sin 2 9* Q B = 0.30tg;& at the 95% CL. With 
the data currently available this is our best estimate of \U e 2\ 2 and is the most accurately 
known MNS matrix element. 

Finally, we have also demonstrated the possibility of probing the solar interior by 8 B 
neutrinos. We have derived a lower bound on the average electron number density over the 
region where the solar 8 B neutrinos are produced which is 50% of the Standard Solar model 
value. 
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